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MODELLING THE BEHAVIOR OF THE HELMSMAN STEERING A SHIP 
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Abstract 

' sapertanker can be considered as a nonlinear system which 
re; ponds very slowly to changes in the rudder position* Moreover 
t h i i type of ships is often unstable in loaded condition, i.e. it 
nas a tendency to start turning either tc the left or to the right. 
The e oroperties make the supertanker very hard to handle. 

I- order to model the helmsman's behavior, a number of tests has 
been performed using a ship maneuvering simulator. The trained 
selects had tc steer a 200,000 tons tanker along a varying course. 
The results obtained from these trials are encouraging* 

A iiscussion of the further research is given. 

I , I n t roduc t ion 

M ost of the investigations concerned with the behavior of the human 
operator as a controller have been executed with reference to pilots 
of aircraft or spacecraft; some work has been done on the control of 
submarines. The human operator as a controller of surface ships, 
however, did not get very much attention until recently, 

Ir the Netherlands in I960 the Institute TNO for Mechanical Construc- 
tion (7NQ-IWEC0) at Delft built a ship maneuvering simulator [ljin 
order to study among other things ship maneuverability, the design 
of nautical instruments and the training of ship crews. In the same 
year Stuurman [ 2 ] executed a series of trials in which he showed that 
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for smell ships the control behavior of the helmsman could very 
well approximated by means of a Describing Function Model* He 
also » ound evidence that for larger ships a nonlinear model prob- 
ably would give a more realistic description of the helmsman's 
behavior. 

In consult with TNO-IWECO it was decided to continue the work of 
Stuurman as a joint activity of the Shipbuilding Laboratory of 
the Delft University of Technology and the Man-Machine Systems 
Group. Special emphasis was to be laid on modelling the helmsman 
of a supertanker with the following goals in mindi To provide data 
on which the maneuvers of this type of shipa under human control 
can be predicted in a number of situations as well aj to enable an 
evaluation of the employment of a human pilot. The results of some 
preliminary experiments have been reported at laat year’s Annual 
Conference on Manual Control [ 3 ], 


The dynamics of a ship depend not only on the nropertiee of the 
ship itself but also on the topology of the surrounding water. 

The motions of a ship in the horizontal plane can be described by 
a set of nonlinear differential equations. Thasa equations describe 
the tranelations of the ship in a direction corresponding to the 
longitudinal axis of the ship and in a direction perpendicular to 
this axis as well as the rotation about a vertical axis through the 
center of gravity. Fig. 1 gives an indication of the variables 
concerned. 


In 1957, Nomoto [ 4 ] showed that for a ship sailing at constant speed 
the relation between rudder angle 5 and rate of turn r can be des- 
cribed by meant of a second order linear differential equation. For 
moat of the smaller shipa this equation gives an adequate description 
of the ship’s behavior in a numbe/ of standard maneuvers. For a super- 
tanker, however, it was found that the behavior was essentially non- 
linear. Based on full scale trials, Bech [ 5 ] proposed to extend 
Nomoto’ s equation with a nonlinear term. This leads to the following 
relation [ft] : 


T t T 2 f(t) + (T 1 +T 2 )f(t)+a 1 r(t)+a 2 [r(t)] 3 « k[t 3 «(t)t«(t)] t 


( 1 ) 
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FIGURE 1: 

The quantities involved in the description of the 
&hip r e maneuvers, 

where r(t) = d*(t)/dt Is the rate of turn, #(t) is the heading 
angle, 5ft) is the rudder angle, and where the quantities a^, a 2 * 

T l* T 2 * T 3 and K ar# constants. These constants are dependent on 
the hydrodynamic properties of the ship, which are for instance 
related to speed, load condition and possible restrictions in the 
surrounding water. 

In this study a particular ship via. a 220,000 tons deadweight 
tanker In loaded condition has been chosen. The principal data of the 
ship were: Length = 310.00 m; Breadth a 47.16 m; Depth = 24.50 mf 
Draft = 18.90 m; Displacement • 238,000 and Froude number s 0*14, 
The constants in Eq. (l) for this ship have been determined by 
Glansdorp [6;7], they are given in Table 1. If a stationary situation 
Is considered, that is, ?(t) s 0, *(t) “ 0, 6(t) s 0, then Eq. (1) 
changes into: 

♦ a 2 r3 s (2) 
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TABLE It 

Constants in the equation describing the relation between 
the rudder angle and the rate of turn for thr supertanker , 
in deep and etill water ; the speed ooneidered ie 7*72 m/eee. 



Fig. 2 represents Eq. (2) ; this static characteristic la given 

for the ship in fully loaded and ballasted condition. The ship Is 
course unstable In loaded conditions it has a natural tendency to 
deviate from the straight course and to start turning either in 
one direction or the other. 



FIGURE 2 : 

Relation between the rudder angle 6 and the rate of turn r 
in the stationary state for the ehip ooneidered in this 
investigation as. found by Glansdorp [?} (a: loaded oondition; 
hi ballasted oondition ), 
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3 The maneuvering simulator 

The simulator consists of a wheelhouse which has the sane appearance 
a3 that of a real sea-going vessel. The fore part of the ship, the 
sea and a coastline are displayed on a screen in front of the wheel- 
house . The total angle of vision of the helmsman is 120°, The 
ir.age of tho fore part of the ship is static, it Is produced by 
two slide projectors which have a fixed position* The coast line 
is generated by means of a point light source and a movable model 
with three degrees of freedom simulating two translations and one 
rotation in the horizontal plane* A blockdiagram of the system. 



FJCUBF 

Block diagram of the TNO simulator. 

Cn an analog computer the dynamics of the ship to be simulated} 
inciuiing th* characteristics of thrust engine and rudler engine, 
can be programed. The computer yields the signals which control 
the «n«'iror,ncn tal display system and also the instruments such as 
compass, rudder position indicator, log, etc. The helmsman has the 
same controls at his disposal for maneuvering as on a real ship 
viz. the wheel, which gives the input to the rudder engine, and 


the telegraph to the engine, which governs the speed of the pro- 
peller* External disturbances simulating the effects of wind, waves 
and currents can also be introduced into the model on the analog 
computer* 

4 The experiments 

In the experiments described here, the simulator has been used as 
a supertanker at full sea in fully loaded condition as well as in 
ballasted condition, moving at a constant spend. The analog com- 
puter was programmed according to £q. (i) based on the constants 
given by Glansdorp as indicated in Table 1. The rudder engine 
was also included in the simulation* its d/remics have been chosen 
according to the Eqs (3) and (i»), 

T 4 6<t) t 6<t) = « d (t)t < 3 > 

|6(t)| < M, (4) 

where is the position of the steering wheel or the desired 

rudder angle* where ? 4 is a time constant of 1 sec and where the 
quantity H is the maximal value of the rotation speed of the rudder 
CO. 045 rad/sec)* 

The subjects wore four trainees of the School of Navigation at 
Rotterdam* The;* were studying for the rank of first or second mate 
after having be-n at sea for several years. All four trainees were 
experienced in steering conventional cargo ships; only one of them 
(subject A) had sailed on tankers up to 90,000 tone dead weight. 
Their task consisted of following a straight course for about half 
an hour, or of following a preprogrammed but unpredictable course 
which changed between +2° and -2° around a certain nominal couree 
for about fourty minutes* This prescribed course acted as a forcing 
function for the man-supertanker system. Because it is not a real- 
istic situation to follow a continuously changing course, a binary 
signal was chosen* The construction of the binary test signal was 
based on the idea that with the experiments to be executed aleo 
linear models of the helmsman* s behavior might be investigated* 
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• Therefore the input signal (desired course) was constructed in 
such a way that about 70 percent of the input's energy was concen- 
trated at only four frequencies. The aeplitude of the binary signal 
was equal to the two degrees earlier mentioned. The four sinusoidal 
components were the third, fifth, eight and thirteenth harmonies 
of a sinusoid with a period of forty minutes; the phase at the 
initial point was chosen at random, 

TABLE 2: Summary of the tee is performed. 


1972 

April 25 


7.DCpm Subject A 

loaded cond. 
course chang. 

Subject B 2 
6 • 00pm loaded cond. 

course keep. 
Subject A > 
loaded cond. 
course keep. 
Subject B A 
9 . C Opr loaded cond, 
course chang. 

Subject A 5 
ball. cond. 
course keep. 
10.00pm Subject B S 
ball. cond. 
course keep. 


1972 

April 26 

Subject C 7 
loaded cond. 
course chang. 

Subject D 9 
loaded cond. 
course chang. 

Subject C * 
loaded cond, 
course ke •»*>. j 
Subject D *~r 
loaded cond. 
course keep. 
Subject C m 
loaded cond. 
course keep. 
Subject 0 
loaded cond.** 
course chang. 


1972 
May 2 

1972 
May 3 

Subject A 13 

Subject D 

ball. cond. 

loaded cond. 

course chang. 

course chang. 

Subject B w 

Subject C 20 

loaded cond. 

loaded cond. 

course chang. 

course chang. 

Subject A ** 

Subjact D M 

loaded cond. 

loaded cond. 

course keep. 

course keep. 

Subject B 19 

Subject C 22 

loaded cond. 

loaded cond. 

course keep. 

course keep. 

Subject A « 

Subject D *> 

loaded cond. 

loaded eond. 

course keep. 

course keep. 

Subject B <• 

Subject C » 

loaded cond. 

leaded cond. 

.curse keep. 

course keep. 


* With disturbances due to waves, 

** With rate of turn indicator. 

earn Dus to a want of time this run could not be performed. 
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In Table 2 a summary of the tests is given. During the experi- 
ments the compass as well as the rudder angle indicator were used. 
Three tests have been executed in which also the rate of turn 
indicator was used. With the exception of two trials no external 
disturbances simulating ship motions originating from wind* waves 
and currents were introduced; these oxternal disturbances^were 
generated by means of a digital computer and consisted of. the sum 
of 23 sinusoids simulating the motions of the fully loaded tanker 
in a following sea with long crested waves. These waves were consid- 
ered to originate from a wine with a force of eight to nine 
Beaufort. 

During the tests the following signals were reoorded co magnetic 
tape i 

e The desired course 
e The eouree of the ship *(t). 
e The desired rudder angle $ d <t). 
e The rudder angle 6(t). 
e The rate of turn r(t) = v(t). 


Modelling the helmsman's behavior 


In the Figs 4 through 6 some examples are given of the time 
histories of the desired course the course of the ship $(t); 

the desL/t! rudder angle the rudder angle 6(t), and finally 

the rate of turn r(t) as recorded during the testB. In all cases 
the records show that the helmsman generates a rudder angle 6(t) 
as output whlc v consists of discrete steps. Hence the records 
indicate that a linear model to describe the helmsman's behavior 
will not fit tba data vary well. To check this presumption, the 
ratio between the energy not located at the four frequencies 
mentioned before and the total energy of the signal consid- 
ered has beet oalculated. The Table 3 shows the results for the 
four subjects, each having performed two rune of forty minutes. 

From the table it can be concluded that the remnant energy, in 
particular of the desired rudder angle is that high in most 

of the runs observed, that it does not seem appropriate to foeue 
the attention on linear descriptions of the helmsman 4 s behavior 
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Survey of the remnant energy if* the signals observed 
for eight tests wi + h four subjects : Loaded conditions 
without rate of t*. '* indicators without disturbances » 


any longer. It suggests that the helmsman bases his decisions on 
when to move and how much to move the wheel on some criterion which 
is, for instance, a function of the heading angle «fr<t); the rate 
of turn r(t); the position of the rudder 6(t), and, of course, the 
desired angle t^(t). 

By the compass the course *(t) is displayed to the helmsman; while 
in the case that no rate of turn indicator is used, it is assumed 
that the helmsman estimates the rate of turn r(t) directly from 
the time history of the course $(t). Based on the records obtained 
by preliminary tests as well as on the records shown in the Figs 
4 * 5 and 6 a model for the helmsman's behavior is proposed 

which spates that if ‘he absolute value of a quantity s(t) which 
is defined as: 

s<!) : [*{») - * d <t>] ♦ CjKt) (5) 

> ~eeds a certain threshold value d^ t the helmsman moves the 
:> tearing wheel into a position 6^ (t) according to Eq. (6): 

4 d *(t) = a* e (t) ♦ bt(t) ♦ c[* e (t)] 3 ♦ d sign * e <t)* (6) 

In this equation the quantity ^(t) = *<t> - * d <t>, and th ® quan- 
tities >(t), $,(t) and f(t) are the values of these variables at 
d 

the moment that the threshold value is exceeded. As a result of 
this action, the quantity )s(t)| wil decrease; if now |s(t)| 
becomes less than a second threshold value d^, it is assumed that 
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Example of the time histories of the desired course tJt)s the 
course of the skip t(t)s the desired rudder angle 6 the rudder 

angle 6<t); the rate of turn r(tj , and the derivative of the rate 
of turn with respect to the time h(t)i subject As without rate of 
turn indioators without disturbances ; loaded - mdition. 
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Example of the t ims histories of the desired course the 

court* of ths ship *(t) ; th* desired rudder angle S.fsJj the rudder 
angle 6( t) ; the rate of tu. t) 9 and the derivative of the rate 
of turn with respect to + m * ,g h(t): Subject B } without rate of 
turn indicator i without u„ . aneee; loaded condition • 


, 649 - 


12 



Example of the time histories of the desired course tJt): the 
course of the ship 4(t)i the desired rudder angle the rudder 

angle &(t); the rate of turn r(t)> and the derivative of tfte rate 
of turn with respect to the time bit): Subject D; without rate of 
turn indicator j without disturbances i loaded condition . 
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the helmsman moves the rudder again Into zero position. 


Parameter Estimation 


The parameters a, b , c, d, c^, ^ and d 2 can be estimated by 
minimizing a quantity defined ast 
i t 2 

J, * f /[*<«)] it. (7) 

0 

where the signal s(t) represents the difference between the human 
operator output $ d <t) and the model output $ d *(t). The minimal 
value of the quantity can be found by partial differentiation 
of this quantity with respect to each of the unknown parameters 
and by setting the result equal to zero. This yields as many equa- 
tions as there are unknown parameters. Btccjse there does not exist 
a simple analytical relation between the input and the output of 
'he model of the helmsman it is not possible to solve the equations 
jwSt-centloned in a simple way. Furthermore, it should bo noted 
that if a model with given parameters should be inserted into the 
control loop instead of the helmsman, this would lead to another 
time Uistory for the quantities $ d (t>, +(t) and r(t). So, in order 
to get an unbiased estimate of the parameters in the human operator 
model, a comparison should be made between the output of the helms- 
man end the output of thu model, where this modal is also part of 
a closed loop system with the ship's model (See Fig. 7 )• The 

parameters can be found by means of a direct search methed mi^i- 
mizxrg t quantity with respect to the model parametar*. With 
the help <ji Eq. (7) the quantity can be defined! 


£ /[*„<*>! 


a 

* * 


The quantity E^ i'dic • v ow well the model output approximates 
the actual output of .ho uelmsman. 
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Application of an error criterion in Buoh a wap that 
unbiased parameters can bo obtained for a model of a 
system in a oloeed loop* 


The correspondence between the time histories of the actual course 
of the ship P(t> steered by the helmsman and those of the ship 
model p*(t) steered by the model of the helmsman's behavior can be 
expressed by the quantity Eg, defined ass 



( 9 ) 


7 Results 

Based on the data of one test, in which the subject had bean 
instructed to ?teer a ship in loaded condition aloag a straight 
course for a period of 25 minutes, the parameters of a series of 
simple models have been estimated according to the para- 
mater estimation method j ns t- mentioned. The first model was a 
linear ones later models were based on the Eqs CS) and <&)• 
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ba it that the nonlinear terms In Eq. (g) ware neglected, The 
results with these sinple rodsls, however* were very poor* end did 
not give eny indication how to modify the mo» * aed la order to 
S*-t o better result. Therefore, the task of tn lmfloar » which 
du ng those teats was to keep a ship in loaded condition on a 
straight course, was changed Into tho task of following a certain 
unpredictable but well defined course. In this case, again tha 
parameters were estimated, and now the roaulto wore much betters 
trey oven led to e direction in which tho structure of tho model 
had to bo changed, finally, the modal as given by the Eqe it) 
ani (6) was obtained. 

The Table *» shows tho values of tho model parameters as well as 
the values o. ths quantities E ^ and calculated from four testa 

TABLE 4: 

Summary of tho parameters of tho mod*, y given by tho 

Eqo (SJ and (6) and tho quontilUe Tj and Fg« 


St'* jset 

A 

B 

C 

4> 

Date 

1972 
Kay 2 

1972 
May 2 

1972 
May 3 

1972 
Hey 3 

Test nr 

15 

14 

20 

21 

a 

3,5 

1,9 

6,2 

2.0 

b [sec.] 

176.9 

155,0 

210.0 

230,5 

c [degr,* 2 ] 

0,0 

0,1 

-0,1 

0,1 

d [degr.] 

1.7 

1,0 

1.0 

1.5 

Cj [ 3SC .] 

20,3 

30,0 

22.5 

30.0 

djtdcgr,] 

•5.0 

0,2 

0.3 

0.2 

d 2 [degv.] 

-5.0 

9,1 

0.2 

0,0 

IjW 

*0,1 

70.9 

43.4 

32.8 

* a M 

8.7 

-- 

i mf 

3,1 


- sweetie the four subjects A, 0, C end 0 (Bee Table 2 >, The 

Figs g .6 g shew some typieel time histories of tbs eet-al 

signals #^(t>, /{*) and i^Ct) compared with the medal outputs 
p"(t>cnd d < |*(t), heme remartee with reapeet to the results obtained 
ean ba made. 

e Although a remarkable difference in tho steering behavior of the 
sub jests exists (See tbs figs e, 6 and * ), tha difference 
of tho parameters among tha dlffarant subjects is rattier small* 
a Xn partieularf cr tha subjects A and D tha model fits the teat 
resulte felrly well, 

e Although the output of the model of tho helmemao's behavior does 
eotalway# resemble tho actual output vary wall tho sourea of tho 
ship generated by tha modal experiment fits extremely well, 
e 2n a few eases, far instance ths test with subject C, ths optl- 
mleati.'n of ths parameters by minlmiaing the quantity d %9 was 
difficult, Ths problem was to find tho abooXota minimum of tho 
quantity J ft in stand of a local minimum 

fu rther research 

Until new, only tha paramatare of tho model have freed estimated 
for four eubjoets steering a fully leaded tanker along a prescribed 
course changing at certain moments from vt degrees to -2 dogreoo 
nnd bash. Also soma teats ware performed with a tanker in ballasted 
condition and with dieturbeneeo (simulating far iaetesoe ths motions 
originating from eeewevec) acting an tha ship, bat at this moment 
parameters have mat yet bean so time ted, s^me of there tarts, in parties- 
iar those with the ship selling In open gee with rather high waves, 
showed that fairly dangerous situations can occur when tha helmsman 
is net wall trained, 

Tha further research is fcouaeed an tha fallowing points which 
may ha of importance in developing a uaafnl modal in order to 
ha able to predict tha behavior of a ship under human control, 
a A g.sdy of tha influence of tha »»ah to ha eueeuted by tha 
tub jests on tha atrueture of tha nodal of tho helmsman's 
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9JQVBB d* Tvpioal tins hi 9 tori ss of ths astual signals p A (t), t(t) and l A <t> oomparsd v< th 
ths mods l or tputs t*(t) and If (t)t Subjoot 0$ loads d oonditionj without pals of 
turn indicator. 






behavior and on the parameters belon ing to the modal. 

• A study of the influence of the ship dynamics on the structure 
of this Model ind on the model parameters* 
o A stuoi of the Influence of the disturbances acting on the ship 
on the structure of this sodsl and on tha nodal parameters, 
a The physiolog leal and j, y ecological lntarpratation of tha 
eodel parameters. 

A large nusber of experiments with tha TRQ ship nanauvaring sim- 
ulator have been planned in the niddla of 1973. Tha anparinanta 
have been eet up in such a way that about tan eubjecte will have 
to steer e large number of ahlpa with different dynanlca. In 
char -crerizlng these ships with reference to their dynamics, three 
groups can be distinguished* l.e.i 

a Course stable shipe with / nor# or leas linear relation between 
rudder angle a and rate of turn r (tea Pig* ICa)* 
a : nirs« stable ships with a nonJinaar relation between rudder 
angle 4 and rats of turn r (dead none) (see fig. 10 b)« 
a c nurse unstable ships with e nonlinear rslstlon bstwssn rudder 
ingle a and rate of turn r (sea rig. 10c). 



F.GURE 10: 

TH a relation between rudder angle 4 and rata of turn r 
for the three olaeeee of ships* 
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A second tspeo t of tha naneererlng dynamise of ohipg is tha tin# 
constant* Each of tha tbros aloaaao of ahps mentioned before will 
be used in oombinstlon with three tine eenotents vio, to, to and 
210 •••• x * this whole area of poosibl# ahlpa la cowered. 

During the toots f ire t -•optioned no diet urbane as eating on the 
ship (wind, waved ) will bo taken irto oaooont. In tbo beginning ef 
1*7* seas additional trUla will is owes w ted in order to cheek tha 
influence of dieturbaneeo noting on the ship on tbo otreetwre of 
the sodsl of the be la scan's behavior and on tha paraaetsrs belonging 
to the eodel. 
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ABSTRACT 

Hunan operator performance in a simulated anti-aircraft 
artillary (AAA) tasfc ia explored as a function of ay aeon and 
environmental paramo tar a* Oparator responses to repeated 
fly-bys are analysed to provide trajectories of both mean 
response and of rssponsa variability. Both the mean and the 
variability of the response ere sheen to chance with chances 
in experimental variables, and comparison of experimental 
data and theoratioal data obtained from the human operator 
model Indicate thee performance chances can be related to 
ehances in tho human* a ability to obtain and prooees infor- 
mation, Xn addition, it appears that the human predicts the 
future course of ths target to a limited extent. Modi fleet ion 
to tho human operator medal to aeoount for this type of be- 
havior are suggested. 
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